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Abstract: A new procedure is presented to develop multi-variable model-based 
parameter estimation (MBPE) model to predict far field intensity of antenna. By 
performing MBPE model development procedure on a single variable at a time, the 
present method requires solution of smaller size matrices. The utility of the present 
method is demonstrated by determining far field intensity due to a dipole antenna over a 
frequency range of 100 - 1000 MHz and elevation angle range of 0-90 degrees. 

I. Introduction 

Recently, there has been considerable interest in the use of the Model Based 
Parameter Estimation (MBPE) technique for improving the efficiency of computational 
methods used in electromagnetic analysis [1], MBPE involves accurately fitting a 
mathematical model (Pade’ rational function with unknown coefficients) to computed or 
measured quantities at a finite number of points. The unknown coefficients of the model 
are obtained from the solution of the resulting set of equations. Miller [2-3] describes in 
detail the theory behind MBPE and gives many examples to which MBPE may be 
applied. In most of the earlier work the Pade’ rational function fitting model is used for 
MBPE in a single frequency variable. In [l], the Pade’ rational function fitting model in 
the frequency domain is extended to include spatial dependence by making the numerator 
and denominator coefficients dependent on the spatial variable. The coefficients in the 
Pade’ rational function are determined by solving a matrix equation with N e xN f 

unknowns, where N g and N f are, respectively, the number of sample points in the 6 
and / dimensions. The technique described in [1], when extended to multi-parameter 
dependence, would require solution of a matrix equation with N g xN j-xN^x.... 
unknowns, where A (<) is the number of sample points and (•) can be replaced by 

6,/,0, In the present communication an alternative method to develop the MBPE 

method for two-parameter dependence is presented. Instead of requiring the solution of a 
matrix equation with N e xN f unknowns, as is done in [1], the present method requires 

solving a matrix equation with N f unknowns N e times to obtain the MBPE model for 

the frequency variable. The MBPE model for the 6 variable is then obtained by solving 
a matrix equation with N e unknown N f times. One of the main advantages of this 

current approach is that the MBPE model for each variable is developed independently 
and hence requires smaller size matrix inversion. 


II. Theory 



Consider a dipole antenna of length L and diameter d 0 placed along the z- 
direction. The radiation pattern of the dipole as a function of frequency / and elevation 
angle 6 can be written as [1] 

\ E \- rif O) _ N o (0) + N 1 (0)f + N 2 (0)f 2 + .. + N n (0)f'' 
ll l + Dl (d)f + D 2 (d)f 2 +..+ D d (0)f d 

where N 0 ,N 1 N n and D v D d I) d arc unknown coefficients and depend only on 6. 

The expression (1) is the required MBPE model with dependence on two parameters. 

The coefficients N 0 , N ] N n and D l ,D d D d can be expressed using 

A nO+ A ,, l 0+\J 2 + 

AM = x + a e+A e 2 + a e* (2) 
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where A n (6) is either N n or D n . The expression (2) is the MBPE models for the 
unknown coefficients appealing in (1). When q = 0, the expressions (2) arc 
identical to expression (7) in [1], The unknown coefficients appealing in (2) can 
be calculated from a set of simultaneous equations foimed by equating (1) to the 
calculated or measured electric field intensity due to the dipole at (n + d + l) 

frequency points and, (p + q + l) points in 6 space. If the procedure outlined in 

[1] is followed, the unknown coefficients are obtained by solving a matrix 
equation with (p + q + l)*(n + d + 1) number of unknowns. In this paper we 
present an alternative approach, which requires solving a matrix equation with 
(n + d + l) or (p + q + l) unknowns. 

Consider the sampled electric field intensity due to the dipole for 6, at (n + d + 1) 
frequency points and using (1), (n + d+l) simultaneous equations are generated. 

The solution of these simultaneous equations are used to solve for N 0 ,N 1 N n , 

D\>Dd D d 

Similarly, from the sampled data for 6f,6 3 , 6 ( 1} at (n + d + l) frequency 

points, coefficients (N 0 (),N 1 ()..N n (),D l (),D 2 ()...D d Q) , at (p + q + l) 6 values 

are calculated. Now considering the coefficients N 0 (6 l ),N Q (6 1 ),....N 0 (6 p+q+l) ) 
and using expression (2), (p + q + l) simultaneous equations are generated and 
solved for (A n0 ,A nl ....A np ,A n(p+1) ,A n{p+2) ,....A, l(p+q) ) . The electric field intensity at 
any 6 and / is obtained by first calculating the coefficients using expression (2), 
and then using expression ( 1). 

III. Numerical Results 

For numerical validation of the procedure described above, a dipole of 
length L = 0.4836 meter and diameter d 0 = 0.0003 meter is considered. Using the 
Method of Moments (MoM), the far field intensity values due to the dipole at 
sample frequencies f s = 0.15, 0.30, 0.45, 0.60, 0.75, 0.90, 0.95 GHz and 
sample elevation angles 6 S =2., 10., 20., 40., 50., 60., 80., 90. degrees are 



calculated. Selecting (n =4 ,d = 2) and ( p = 5, q = 2) various coefficients for 
MB PE arc determined. The far field intensity for any □ (0<6 <90 degrees) 
over the frequency band 0 < / < 1 .0(1 Hz. now is obtained from the knowledge of 
(N 00 ,N 01 ,....N 0(p+q) ), (n + d + 1) ..., etc., coefficients. Figures 1(a)- 1(e) 
show far field intensities calculated at 6 = 10, 30,50; 70ctnd90 degrees over the 
frequency band using the MBPE model. Also plotted in Figures 1(a)- 1(e) arc the 
far field intensities calculated using the MoM. As can be seen from the figures, 
the results from the MBPE method presented in this paper arc indistinguishable 
from the MOM calculations. In Figure 2, the far field intensity for this dipole as 
computed by the MBPE method presented in this paper is shown as a function of 
both frequency and theta. 

IV. Conclusion 

A new method has been introduced in this communication to develop 
multi-variable MBPE model to predict an antenna radiation pattern as a function 
of frequency / and spatial angle 6 . One of the main advantages of the new 
method is that the order of the matrix to be solved to determine MBPE parameters 
is much smaller than the earlier approach [1], 
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Figure 1. Far field intensity as a function of frequency for a dipole with length L = 0.4836 meters and 
diameter d 0 = 0.0003 meters as calculated by MOM and MBPE at theta = 10, 30, 50, 70, and 90 
degrees, with the seven fitting frequencies shown as points. 







